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Abstract 

We present an improved method of reconstructing New Physics (NP) masses 
from invariant mass endpoints. While the traditional method focuses on a sin- 
gle NP decay, our method considers the decays of two or more NP particles 
{ABC.) in a grander decay chain: anything — > ABC... ... — > jets + leptons. 

Though the center-of-mass energy Ecm of 'anything' varies unpredictably at 
a hadron collider, a sample of many events nonetheless expresses features of 
threshold production Ecm = mA + mB + .... invariant masses constructed from 

the final jet and lepton momenta are correlated in a way that makes their thresh- 
old endpoints visually obvious in a scatterplot. We illustrate this technique for 
the production of two neutralinos in the MSSM: anything — > XiX^ {i-iJ = 
2,3,4) which subsequently decay via on- or off-shell sleptons to four leptons. 
Assuming the relevant SUSY spectrum is below 1 TeV and squarks/gluinos 
eventually decay to neutralinos, our MC study shows that one low-luminosity 
year at the LHC (10 — 30/6~^) can quantitatively determine on- versus off- 
shell decays and find the relevant neutralino and slepton masses to less than 10 
percent. 
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1 Introduction 



At a hadron collider such as the Tevatron at Fermilab or the Large Hadron Col- 
lider (LHC) at CERN, searches for New Physics (NP) states beyond the Standard 
Model (SM) must take into account the fact that partonic center of mass energies at 
these machines are not tunable (as they will be at a much anticipated e^e~ linear 
collider [1]) but vary continuously in principle from zero to the combined hadronic en- 
ergies of approximately 2 TeV and 14 TeV, respectively. Moreover, many NP models 
predict the production of a long-lived particle that is likely to escape the detectors, 
carrying away missing energy. The traditional method of seeing a NP mass as a sharp 
resonance in a cross section is therefore not applicable, and we must consider other 
approaches to precision measurement of NP masses crucial for testing properties of 
an underlying fundamental theory. 

One well-studied avenue is to construct invariant mass distributions of final jet 
or leptonic momenta in exclusive channels and study their endpoints. Even if some 
NP particles carry away missing energy in each event, endpoints of said distributions 
can be measured and matched to analytical functions of NP masses [2]. There are 
several caveats to this method however. First, the exclusive channel under study 
must somehow be identified or assumed. Second, backgrounds must not interfere 
with endpoint measurement. Third, there may be some model-dependence in the 
method of fitting the endpoint on a 1-dimensional histogram. The first caveat is the 
most difficult to deal with, especially in a model such as the Minimal Supersymmetric 
SM (MSSM), where the gluino and squarks decay via literally hundreds of possible 
decay chains(" cascades"). Studies of cascades which have enough endpoints to solve 
for MSSM masses, e.g. g qq ^ X2II ~^ l^l^qq l^l^qqXi, usually just focus on 
a few 'benchmark scenarios' jSl S] which, however, may not be what Nature chooses. 
As for eliminating SM backgrounds, requiring a suitable number of hard jets and 
isolated leptons may suffice; NP backgrounds are more challenging and, if these can 
be reduced, typically also inflict damage in the region of the desired endpoint where 
rates are already This then brings up the third problem of how to fit the 

endpoint. Linear or Gaussian fits are most convenient, but very detailed study of 
cuts and detector effects are required to understand their general accuracy [5J. 

In this work we wish to show that there are two important features of NP particle 
production which, when used together, can greatly boost the efficacy of the endpoint 
method. First, if NP particles carry a new conserved charge, such as R-parity in 
the MSSM, they will be multiply-produced. We may, for example, consider inclusive 
decay chains of the form X AB jets -|- leptons, where A and B are NP states 
and X is any system of particles with a sufficiently large total invariant mass: mx > 
niA + TTiB (the case where mx = mA + ms is especially important and we designate 
this 'threshold production.'). Since X usually includes at least one multi-particle 
system, mx really is continuously-valued from the threshold value all the way up 
to the machine energy. Second, depending on mx and the exact way in which AB 

^One notable exception is the decay A ^ B ^ C l^P, where A,B,C are NP particles, e.g. 
in the MSSM ^7 Xi- Here the dilepton invariant mass distribution rises up to its endpoint and is 
therefore robust in the presence of diffuse backgrounds. 



decay to the specified endstate, invariant masses constructed from tlie final jet and 
lepton momenta attain endpoint values for certain kinematic configurations only. 
At threshold production, in particular, one special configuration will simultaneously 
maximize several invariant masses; collecting a large number of threshold decays, 
a 2-d or 3-d scatter plot of these invariants would exhibit a clustering around this 
'threshold point.' Yet threshold production is clearly only an infinitesimal possibility 
and superimposing events for 'above-threshold' production (m^ > itia + trb) would 
probably hide the threshold point (hereafter called a 'hidden threshold'). Contrary 
to this intuition, however, we find that, for some invariant mass combinations, the 
hidden threshold is highly visible, being in fact fortified by above-threshold events. 
This allows us to directly measur^ values of threshold invariant mass endpoints 
which, as usual, can be used to constrain NP masses. Moreover, backgrounds should 
not be a problem since these have different invariant-mass correlations. 

The precise mechanism of this Hidden Threshold (HT) technique is best illus- 
trated by example, which in this work we take to be 

X - - i^rx^) i^,J = 2,3,4) (1) 

i.e. neutralino pair-production (via any parent channel) and decay to leptons {i = 
e,/i). Section [2] explains the basic HT theory in this case, followed in Section [3] by 
application to Monte Carlo (MC) generated data simulating a low-luminosity (10 — 
30 fb"^) run at the LHC for two different MSSM parameter points: one with on- 
shell slepton decays, and another with off-shell decays. Finally, Section |4] summarizes 
these results and suggests many avenues for further application. 

2 The Hidden Threshold Technique 

2.1 Base Case: Pure Higgs Decay 

Let us explain the HT technique in three steps, the first of which will be the case 
where X in ([1]) is a single particle which therefore has a fixed mass, say a heavy Higgs 
boson {H or A, hereafter designated 'Higgs'): 

pp ^ H/A ^ x° ^ e^e^fi^fi^ ^ e+e-/i+/i-x?X? (2) 

For now, let us take the intermediate sleptons to be on-shell (see Sec. l3.2l below for the 
case of off-shell decays). We studied this decay chain in [5], and for the convenience 
of the reader we briefly recapitulate. From the four observable final lepton momenta 
in ([2]), one can define seven independent invariant mass combinations. These we 



■'There is also no need to 'fit' these endpoints, as we will see later. 



defined as follows (leptons labelled as 1, 1', 2, 2', see Fig. [T]): 



Ml = {PI+PV +P2+P2'f (3) 

m\i2i = {{pi + Vv - V2 - V2'Y + {pi + V2' - V2 - Vvf (4) 

+ {Pl +P2- Pi' - P2')^}/3 

Mm = {{Pi + Pr + P2 - P2'Y + {Pi + Pi' + P2' - P2Y (5) 

+ {P1 +P2 +P2' -Pl'Y + {P2 +P2' +Pl' -pi)^}/4 
^2/ = {{Pl + PV - P2)^ + {Pl + Pi' - P2')^ + {Pl + P2 - P2')^ (6) 
+ {P2 + P2' - Pl')^ + (Pl - Pi' + Ps)"^ + {Pl - Pi' + P2')^ 
+ {Pl -P2+ P2'Y + (P2 - P2' + Pi')'' + {PV +P2- PlY 
+ {Pl' +P2' -Pl)^ + iP2+P2' -Pl)'^ + iP2' +Pl' - J52)''}/12 
Mil = {{pi + Pi' + P2y + (pi + Pi' + P2'Y (7) 

+ {Pl +P2+ P2'Y + {P2 + P2' + Pl')^}/4 
Ml ^ {(p,+p^,)4 + (p^+p2')' + (Pl+P2)' (8) 
+ {P2 + P2'Y + {P2 + Pl'Y + (Pl' + P2'Y}/Q 
04 = PlPvP2P2'(^IJ^upa (9) 

Adding the usual dilepton invariants Mee and M^^ to this list, we studied all 
distributions (1-d histograms) and derived analytic formulae for their associated 
endpoints (see ^ for exact expressions, too lengthy to reproduce here). In MC 
simulations of LHC data (assuming a luminosity of L = 300 fb~^) at several MSSM 
points where both neutralinos and both sleptons were degenerate, endpoint preci- 
sions were not high enough to give better than 30% determination of the unknown 
masses niHi ^x°' ^^"^ "^x? (hereafter we abbreviate irts = mj, irti = m^o). This 
rather unexpectedly poor resolution from a seemingly over-constrained system(sever0 
constraints on four unknowns) owes to the highly non-linear form of the endpoint 
formulae, which generically yield a discrete set of solutions — adding in smearing 
and detector effects melds this discrete set into the large continuous range quoted. 
However, if one of the masses were already known to 5% accuracy, the degeneracy of 
the solutions breaks and the other three masses are likewise well-determined. Though 
our choice of invariants in ([3])- ([9]) is somewhat more complicated (yet more system- 
atic) than is customary, the above is nonetheless a standard illustration of the usual 
invariant mass endpoint method. 

2.2 Threshold Higgs 

Next, we consider what happens if the two neutralinos in ([2]) are produced at thresh- 
olcl^, that is when uih = rui + rrij exactly. Now, it turns out, five of the invariants 



^We did not derive an endpoint formula for a4; however, we could numerically see it did not make 
a difference. 

^This situation was far from realized in 8J, where the Higgs mass was set fairly high, rriH ~ 
400 - 600 GeV. 



defined above (M4;, M2i2h ^nu M^i, and Mn), in addition to M^e and M^^, are si- 
multaneously maximal when the kinematic configuration shown in Fig. [T] is realized. 



Figure 1: Kinematic configuration which simultaneously maximizes M41, M2121, M121, M^i, 
Mil, O'f^d (Mee X M^^) at threshold (mn = mi + ruj). 

With a suitably large sample of events, we need only plot one invariant against 
another to get its maximal value, as shown in Figure [2l Note that the endpoints tend 
to lie at a fairly triangular apex, therefore well-approximated as the intersection of 
two tangents to the 'correlation shape'. Moreover, this correlation retains its shape 
in the presence of backgrounds, since these latter do not have the correct correlations 
among the various invariants and would tend to form a diffuse halo around the more 
concentrated signal shape. Compared to the traditional one- dimensional histogram 
approach, where backgrounds are harder to subtract and endpoints must be fit with 
a more arbitrary function, the reader can begin to appreciate that a technique using 
correlated invariants is much more powerful. Though this case of threshold produc- 
tion is not particularly likely in the MSSM, it prepares us for the next step. 

2.3 Continuous Superposition of Higgses 

Consider now the case where not a single Higgs decays in ([2]), but rather a continuum 
of Higgs with masses in the range mj + mj < mn < 00. Fig. [3] shows what to expect 
for several strategically-choseiij choices of invariants: threshold decays are plotted in 
gray, while all above-threshold decays are plotted in blacl^lll. The key observation to 
make here is that the threshold points P(Q) in Fig. [3^(b) are not at all obscured when 
above-threshold events are superimposed; in fact, these latter serve to graphically 
reinforce the threshold points, located on the envelope of the collective shape. 

It should be clear, however, that without our color-coding in Fig. [3](a,b), it is not 
possible to tell the exact position of P and cj^. Nevertheless, P and Q share a tight 
relationship: namely, they must uniquely identify the three endpoints M^"-^, M^i , 
and (M™"^) x (MJ^^^), this last of which is obtainable from a wedgebox plot |lUj . 
i.e. a plot of Mee vs. M^^. We can therefore find the precise locations of P and Q 

^For the specific decay topology ([T]), it turns out that only certain correlations between invariants 
axe easily analyzed, these bein_g_M4M^s. MaiMgeM^^, M3/ vs. M^iMeeM^^, M31 vs. M41M2121, 
M31 vs. MaiMi2u and Mn vs. M2mMm. 

^Threshold decays obviously contribute only infinitesimally to the total shape, but for the sake 
of seeing how these are distributed compared to above-threshold decays, we plotted 10'^ of these on 
top of 10^ above-threshold events. 

^There is a sort of kink near P and Q, but this turns out to be related to M]"'^f: only. 
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Figure 2: Various correlations of four-lepton invariant masses at threshold 

(ruH = rui + nrij). This simulation is purely from relativistic kinematics with 
(mi, rrii, mj, mg) = (1, 3, 4, 1.1) (units arbitrary). Dashed lines mark the maximum 
value of each invariant. 



with a graphical device: superimpose a plot of M41 vs. M31 ^jl^mTJ^y.\^jJlilL) with one of 

(M™^=^)x(M™'^==) (i-e. superimpose Fig. [3^ and Fig. [3b with swapped and 

rescaled axes). This gives rise to a characteristic 'cat-eye' shape (Figure [3b), where 
the upper corner of the eye pinpoints the correct threshold extrema of M41 and M3/. 

There is another very useful correlation here, namely M^i vs. M4iM2m (see pre- 
vious footnote for alternatives). Analogous to the situation in Fig. [3^ and Fig. [3)d, 
threshold values of M3; and M^°'^M2i2i for this X — Xi'Xj decay again lie on the 
envelope of the correlation shape at point 'X' in Fig. [311; what is more interesting, 
every point on the envelope above X (say, 'Y' in Fig. [311) corresponds to another set 

of Mgf and M^"-"" M^l^i for a decay X ^ C^, where C and ^ are some hypothetical 
set of heavier particles (m^,^ > rrii^j) which follow the same decay chain through an 
on-shell slepton to leptons and the LSP. In other words, if ( and ^ really were pro- 
duced, they would have decay kinematics leading precisely to the endpoints identified 



at Y (we will call this a 'pseudo-threshold'). Every such pseudo-threshold therefore 

yields a set of endpoints Mg; and M^"^ ^2«2« which can be used to constrain mi, 
TTT-s, and (where without loss of generality we can take C = 0) ^^is latter being 
merely an auxiliary parameter at each such point. What we are effectively saying 
here is that the shape of the envelope above X depends on mi and m^ only, and 
sufficiently precise measurement of this envelope can constrain these parameters. 




Figure 3: Correlations of various invariant masses in a superposition of many Higgs decays 
with a continuum of masses in the range mi + mj < mn < oo ( again from relativistic 
kinematics only, with (mi, mi = mj, mg) = (1, 2, 1.5) in arbitrary units). Points 
P and Q in (a) and (h) identify threshold endpoints; these can he found precisely with a 
'cat-eye' plot in (c), where we have abbreviated [3 = (M™"^') x (M^"^). In (d) we show 
another useful correlation of invariants which not only gives threshold endpoints at X, but 
also 'pseudo-threshold' endpoints at any point Y higher up on the envelope. 

Now there is really no difference between this situation with a continuously- 
massed Higgs and that of the decay ([1]), or more specifically 

X ^ X' + x°(- e^e^ -> e+e-x?) /^^/^^ - /i+Z^-X?) (10) 

where X could be pp, qq' , qg , gg, X2X2 ^tc, or any twc@ SUSY particles with a 
continuous center of mass energy greater than or equal to mj -|- ruj, while X' are any 



The present work assumes R-parity is exact, though the general technique does not require this. 



collection of particles, e.g. hadronic jets, that do not confuse the 4-lepton signal. 
As long as the neutralino-pair subchains are intact, the mother chain is irrelevant. 
Note also we do not require the lightest neutralino to be stable, as long as its decay 
products do not include leptons. As we shall see in the next section this technique 
works very well in MC simulations with a minimal requirement of four isolated leptons 
plus missing energy. 

3 Examples of Applications 

In this section we demonstrate the robustness of our technique at two different MSSM 
parameter points with varying decay topologies. Our MC setup uses HERWIG 6.5 
and private codes as in our previous publications [Bl El [101 reader is referred 

to these for details. 
We generate events 



for 10—30 fb~^{a low-luminosity year at the LHC) and only pass those which decay to 
a hard and isolated e~^e~ and pair with sufficient missing energy(i^7^ > 20 GeV). 
This eliminates SM backgrounds aside from pp — > Z*Z, though this can be modeled 
and subtracted (by a Z-veto if necessary, which happens to work perfectly at the pa- 
rameter points we will consider). SUSY backgrounds fall into two categories: those 
which produce exactly four leptons and those which produce more than this (presum- 
ably losing some to isolation cuts, detector effects, etc.). The first category includes 
slepton or chargino '3+1' decays such as xf (— > l^l'^l'"^ ui/'u' Xi) Xii^ ^"^ ^" Xi) 
where one decay goes to three leptons and the other to one. These could always be 
substantially reduced via fiavor subtraction, but this is not necessary since they have 
totally different kinematics from a '2+2' decay (i.e. (ITO!) ) and should not obscure 
the envelopes in plots such as Fig. [31 Moreover, such decays would exhibit their own 
characteristic correlation shapes in addition to the trilepton invariant mass edge. In 
the worst case, then, four-lepton SUSY backgrounds should be considered an 'en- 
riched signal.' The second category of SUSY backgrounds includes events such as ff 
where each stau decays r ^ l^l'^l'^ ^^'Xi two of the total of six leptons 

are somehow lost. Yet such events, if not made utterly small by leptonic branching 
fractions, would introduce a more-or-less diffuse halo in our correlation plots, and in 
particular should not confuse the identification of threshold points. 



The most basic illustration of the HT technique for the decay ([TOj) is when only one 
particular (i, j)-combination has a significant branching ratio; if, in addition, i = j, 
threshold endpoint formulae simplify dramatically, e.g. 



PP {<?, 9, Xt2^ X?,2,3,4} + {Q^ 9, Xt2^ X?,2,3,4} 



(11) 



3.1 On- Shell Box 



Ml 



max 



m'j{2m'jm1 — ml) — 



(12) 




Table 1: Relevant masses at the On-Shell Point and Off-Shell Point (all masses in GeV). 









xH 


xl 




xt 


X2 


On-Shell 
Off-Shell 


117.1 
79.6 


197.5 
131.5 


257.2 
160.6 


317.7 
249.5 


134.7 
148.0 


193.3 
115.6 


317.0 
249.0 



A X2X2 pair, for example, is often the chief product of colored cascades at mSUGRA 
points such as SPSla [3]. Though we have checked our technique works at SPSla, 
we choose another point which has higher rates and therefore better demonstrates 
the agreement with theory presented in the last section: 

On-Shell Point 

/i = 250 GeV M2 = 250 GeV Mi = 125 GeV 
tan /? = 10 Mr,[e,-^)^ = 250 GeV M(e,^)^ = 130 GeV 
tuh = 700 GeV Mg ^ 400 GeV Mg ^ 500 GeV 

Here we have raised the left-handed slepton and stau soft mass inputs above those 
of the right-handed selectron and smuon (note the physical masses differ slightly 
from these, see Table [1]) to suppress sneutrino exchange and stau modes which may 
reduce the magnitude of (but not character of) the signal. The heavy Higgs mass is 
set rather high to remove it from the analysis, but setting it lower would in fact be 
beneficial since Higgs contribute to the signal process X —>■ x^x^. 

With a luminosity of 30 fb~^ the wedgebox plot in Figure H^i is a very dense and 
symmetric box structure, which strongly suggests degenerate neutralinos decaying 
via degenerate slepton^lj. Since the dilepton mass distribution(Figure Hb) is quite 
triangular, we could also guess that the sleptons are on-shell — in the next subsection 
we will show how to establish this quantitatively. 

Thus, being certain of having a fairly homogenous sample of X — X2X2 
cays, where each neutralino decays as X2 ^^^^ ^^^~X2 with a dilepton edge 
Mi+i- ~ 75 GeV, we can directly proceed to construct a cat-eye plot in Fig. Hl^ to 
determine threshold values of M^"-^ and M^i . Note that even at this low luminos- 
ity, the threshold point (marked as a large black dot in Fig. St) can be located to 
a few GeV precision: we measure (M^"^, M3; ) = (150 ± 5, 107 ± 5) GeV. This is 
quite close to the expected values of (153.9, 112.4) GeV found by plugging in masses 
from Table [1] into formulae in the Appendix. The HT method is working. 

Having determined M^i from the cat-eye plot, we proceed to measure a number 
of pseudo-threshold points in a plot of M31 vs. M^iM^m for M^i > M^i (see Fig. Hi) . 
Table [2] compares some of our measurements to expected values over a liberal range 
of M^i, where the reader can verify good agreement. Using a large number n of 

^"^Strictly speaking this could also be a chargino decaying via a sneutrino, i.e. xt ~^ ^ 
I'^l'^xt- Other signs such as the number of 6-lepton events (after xt ~^ l^'^Xi) rnight resolve this 
ambiguity. 
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Figure 4: Various plots for the On-Shell Point with 30 luminosity. The wedgebox plot 
(a) is mostly box-like with a triangular dilepton mass distribution (b) typifying an on-shell 
decay. The threshold values of M41 and M^i for the X2X2 resonance ((3 = (75 GeV)'^) in the 
cat-eye plot (c) are easily identified at the 'corner of the eye' (the black dot). We proceed 
to measure several points along the envelope of the shape in (d) for M^i > M^i ; these can 
be used to constrain the slepton and LSP masses. 



measurements (say n = 30) we search (mi, m^, mcijUir^, ...mr„)-space for the best fit 
to the corresponding (Mg^ , M^°'^M2i2i )i=i..n, using analytical expressions in the 
Appendix. We find rui = 116 ± lOGeV with the slepton mass shifted by a near 
constant, rris = (mi + 17) ± 2 GeV. From these and the dilepton edge M^+i- = 75 ± 
1 GeV, we can likewise determine m2 — 210±12 GeV. With a higher luminosity 
sample (say, 300 fb"^) error-bars could certainly be reduced by a factor of several; 



Table 2: Measurements of M41M2121 for various M31 at the On-Shell Point (viz. Fig. 
and comparison to theoretical values (boldface) assuming these correspond to pseudo- 
thresholds, i.e. plugging slepton and LSP masses from Table {1\ and (in place of mj) 
into formulae in the Appendix (all masses in GeV). 
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318.0 
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363 ±2 
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235.2 
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269.0 
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875 ±5 
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313.0 


300 


1260 ± 5 


1250.6 


595.5 



in this case the masses (mi,ms,m2) can ultimately be found to about 2% accuracy. 
3.2 Off-Shell Box 

Decays through off-shell sleptons are simpler for two reasons. First, mathematically 
speaking, invariant mass edges do not depend on slepton masses; in the case of a box 
topology these can therefore only depend on mi and mj. Since the dilepton mass 
edge is equal to the difference of these, {nij — mi), we only need one other function 
of these masses to determine them. 

Secondly, the physical degrees of freedom are easier to analyze: for threshold 
decays in particular, each neutralino x'j is at rest in the center of mass frame and 
leptonic invariants are maximized/minimized when the same flavor leptons are emit- 
ted antiparallel/parallel. The product (Mge) x (M^^) is therefore maximal when the 
electron(muon) is antiparallel to the positron(anti-muon), and is equal to (m^- — mi)^. 
For this kinematical configuration, however, all the other invariants M4;, M2121, M121, 
M31, and Mil are equal(up to a constant factor) to (m^ — mi) and therefore do not 
provide independent information. 

The trick is to consider configurations where the electron and positron are an- 
tiparallel while the muon and anti-muon are parallel (or vice versa). In this case Mee 
or M^^ is maximal and the other invariants (excepting M121) unconditionally attain 
minima when the e+e~ pair is emitted at right angles to the /i"*"/!" pair: 

-rr-jmin 

^41,2121,31,11 

for specific values of the parameters a, /3, 7, and ^ listed in Table [31 Thus, any 
one of these in conjunction with the dilepton mass edge uniquely determines mj 
and mi. Geometrically, these endpoints are found at the intersection of the line 



{nij — nil) 



(13) 



Mi+i- = {rrij — nil) with the near edges of the correlation shapes on a plot of Mi+i- 
versus M4;, M2121, M^i, or Mu. These intersections can be found with high precision 
and then checked against eachother for consistency. 

To see the efficacy of this method, consider the following parameter point (see 
Table [1] again for the physical spectrum) : 

Off-Shell Point 

/i = 150 GeV M2 = 200 GeV Mi = 100 GeV 
tan ^ = 10 M^,(e,/i)^ = 300 GeV M^^,^)^ = 135 GeV 
ruH = 700 GeV Mg ^ 300 GeV Mg ^ 350 GeV 




Figure 5: MC Simulation at the Off-Shell Point for 10 fb~^ luminosity. The dilepton 
invariant mass distribution in (a) indicates a maximum of about M™^5 ~ 52 GeV. In (b) 
we show where this maximum is coincident with the minimum value of M2121 CLt the dotted 
line, which is to be matched against expression U3\) . Corresponding plots for the other 
invariants look very similar. 

Though the 10 fb~^ dilepton mass distribution in Fig. [5K looks very similar to 
that of the On-Shell Point, the dilepton edge at M^I^T: ~ 52 GeV is in fact due to 
three-body decays via off-shell sleptons. We might guess this from the vaguely less- 
than-triangular shape, but now we have a better way. We simply plot Mi+i- versus 
each of the invariants M^i, M2i2h M31, and Mu (e.g. Fig. EJ^b)) and match against 
the formulae in ( |T3i) . Table [3] shows that these give mutually consistent values of 
mi and ruj, which at high luminosity may be determined to several percent or so. 
Had we constructed the same plots at our On-Shell Point, we would have obtained 
a set of endpoints which give inconsistent, even nonsensical values for mi j (e.g. 
negative LSP mass!). Though the dilepton distributions of Fig. |Dd and Fig. give 
qualitative hints of whether sleptons are on- or off-shell, here we seem to have the 



Table 3: Offshell Decay Parameters in il3\) and Fits to MC data at the Off-Shell Point for 
the LSP mass (mi), assuming the dilepton edges are exact and 0.2 GeV precision on other 
endpoints. Also shown are attempted fits to data at the On-Shell Point which, even with 
larger(zizl GeV) errors, give self-inconsistent results. 
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first definitively quantitative method of determining this which succeeds even with 
modest (several GeV) endpoint precision. 

4 Conclusions 

Let us now summarize the HT technique and how we applied it in this work: 

1. DEFINE the decay chain X ^ X' + ABC... —>■ rij jets + ni leptons, where X is 
anything and X' is exclusive of jets and leptons. In our case this was X — x^Xj 
where each neutralino then decayed via a slepton to a pair of leptons and the 
LSP. 

2. DERIVE analytic expressions for all jet and lepton invariant mass endpoints as 
functions of NP masses, and find kinematic configurations for which two or more 
of these are extremal at threshold production. In our case of four leptons(ni = 4 
), there are in principle 7 independent invariants; of these, we found (Mee) x 
(M^^), M4;, M2121, M121, M31, and Mil were simultaneously maximal for on-shell 
decays, wheareas for off-shell decays Mee(M^^) was maximal(minimal) where 
M41, M2121, M^i, and Mu were minimal. 

3. DISPLAY correlations of the above invariants in a scatterplot which makes their 
threshold endpoints visually obvious. Use these to solve for NP masses. For 
off-shell slepton decays it sufficed to plot each invariant versus Mi+i- and note 
the intersection with the line Mi+i- = MJ+p . For on-shell decays we found it 
useful to make a cat-eye plot to find the threshold value ofl^i'^^ , followed by 
measuring several pseudo-thresholds in a plot of M 2,1 vs. MiiM2i2i- Matching to 
analytical expressions constrains {mi,m2,ms) to within a few percent of their 
nominal values. 

Let us add a few remarks on these three steps. The first step depends of course on 
the specific NP model under consideration, but once a decay chain has been selected, 
deriving analytical formulae for endpoints in the second step is not difficult: the 



sought-after kinematic configuration is usually where each jet or lepton is emitted 
at a polar angle of 6' = or 6' = tt or perhaps 9 = ti/2 (for minima) in the frame 
of the decaying parent particle, backwards Lorentz-boosting to the center-of-mass 
frame (see Appendix of ^ for an example). Choosing which invariants to plot in 
the third step is a matter of trial-and-error — kinematic simulation (as in Fig. [2] and 
Fig. [3]) can be used to see which correlations have a visible threshold point when 
above-threshold decays are superimposed. 

To demonstrate the general applicability of this technique to NP, we present 
below a partial spectrum of examples. Each example in itself entails numerous vari- 
ations(e.g. on-shell intermediate states could likewise be taken off-shell): 

• Neutralino Decay via higgs: X2X2{~^ x\h{—^ bb)). This nj = 4 decay is a 
potential competitor to the sleptonic modes considered in this paper. 

• Sneutrino Pair Production: uu[-^ I'^Xii.'^ Xi^^ (^ ^'^^'))) Here n/ = 4 as 
for neutralino pair production, but with totally different kinematics: one can 
analogously define M4/, M2121, Mi^i, etc., though threshold extrema of these will 
have a different correlation. 

• Chargino Pair Production: xfx^i^ /^z/(— I'^Xii^ Puxi))) This is a n/ = 6 
case which means there are as many as 45 relativistic invariant 

• Slepton Pair Production l^X^i^ ^"^^^Xi))- This could either be = 6 
or = 4 depending on how the other slepton decays. 

• Squark Pair Production gg(— qX2{~^ I ^^Xi)))- ^ = 2, = 4 decay 
with 45 invariants to analyze. 

• Gluino Pair Production gg{-^ l^i^ Q'X^l"* I'l'Xi)) This rij = 8 decay has an 
astounding 588 invariants. 

• Nonstandard Higgs decays h — > Xy(e.g. [TT]) 

• Exotic Vectors (from, e.g. Little Higgs [12| ) 

• Kaluza-Klein Pair Production (e.g. [13]). 

• Exotica decays to top-pairs pp ^ X ^ ti (e.g. [I4j) 

• Low-Scale Technicolor (e.g. W^ttt l^vhb [15j). 

The hidden threshold method provides another way to see signatures of these 
NP models and determine any unknown masses involved. Results can be made even 
stronger in combination with inclusive techniques and complementary methods of 
mass determination (in SUSY for example, see [6l [HI [IZj). We hope that research 
along these lines will allow an earlier discovery of NP at hadron colliders. 

^^The number of invariants is computed as the number of pairwise contractions among the ob- 
servable momenta, as well as with powers of e^"''"' 
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Appendix 

Threshold maxima for on-shell decays are readily computed via the methodology of 
[8]. For the case i = j we just substitute rriA = '2mj in the formulae of that work for 
the [-1- H ] configurations, obtaining the following: 

jyjmax ^ m|(2m|m2-m2)-mf 
41 nijm'j 

]^2l2r = 31/4^ ,^2 (3mfm^ — 4m5m^m^(2m| + m^) — imfm'jm^iSm^ — 12mjmg + 6mjm'^ + 
m^) -|-m^ (16m^ — 32m^m^ + 24m — Sm'jm^ + 3m|) + 6m|(4m^m^ — 4m^m^ + 3m^m^))^^* 

Mj^f^ = Qi/iln rrii (3?^i^^ — 2mfm^ml{5m'j + ml) + Qm\m^m^{?>m'j — m'^rn'^ + mf) — 
2mfmjmf (8mj — Gm^m^-I-Sm^mf +mf )-|-mf (8m^ — IGm^m^-t-lSm^m^ — 10m|mf-|-3m^))^/^ 



max (m^m^— 2mjm^m^+2m^m*— 2m|mf +mf )^ ™j ~ +™s 

^31 = 2l/4m,.mi 



max A /m|m^— 2m^mjm^+2m^m|— 2m^m^H-m^ 

"iwr _V ^ 

11 ~ 24i/4m^m2 



